Tight-binding Hamiltonians for Sr filled ruthenates: application to the gap anisotropy 

and Hall coefficient in Sr 2 Ru0 4 
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Accurate orthogonal tight-binding Hamiltonians are constructed for ferromagnetic SrRuC>3 and 
the layered perovskite superconductor, Sr2Ru04 by fitting to all-electron full-potential local density 
band structures obtained by the linearized augmented planewave method. These Hamiltonians al- 
low the band structure to be computed on very fine meshes in the Brillouin zone at low cost, and 
additionally have analytic band velocities, while retaining the accuracy of the full-potential elec- 
tronic structure calculations. This greatly facilitates calculation of transport and superconducting 
parameters related to the fermiology. These features are exploited to calculate the Hall coefficient 
■ and vortex lattice geometry for Sr2Ru04 with fine integration meshes. We find the lower limit for 

the interband order parameter anisotropy to be compatible with the observed square geometry. We 
also find that the sign reversal of the Hall coefficient can be explained in a conventional way if the 
' bands are shifted by a few mRy so as to match the experimental de Haas-van Alphen areas exactly, 

and the temperature dependence of the relaxation time is strongly dependent on the band character. 

INTRODUCTION 

The Sr-Ru-0 system initially attracted interest because of technological applications as perovskite substrates. In 
the last few years, however, the unique electronic properties of ruthenates made them objects of interest from the point 
of view of fundamental solid state physics. The root compound for this system is the pseudocubic perovskite SrRuC>3. 
This material has an orthorhombic Pbnm (GdFeOa) structure, and is a strong ferromagnet. It is characterized 
by unusual transport properties (so-called "bad metal" behavior), unusually strong coupling between the spin and 
£j ■ charge degrees of freedom, and substantial involvement of oxygen in its magnetic properties. Sr2Ru04 occurs in a 
body centered tetragonal K 2 NiF4Structure and is superconducting at a temperature of 1.5 K. It is generally accepted 
now that the superconductivity of Si~2Ru04 is unconventional and most likely triplet. A number of researchers 
expressed an opinion that ruthenates, or at least some of them, may be strongly correlated systems close to Mott- 
. Hubbard transition. On the other hand, unlike common strongly correlated materials (3d oxides, high-T c cuprates), 
the physical properties of ruthenates are well described by conventional density functional theory, and their unusual 
transport properties are likely to be due to strong magnetic interactions and magnetic scattering. n 

Density functional electronic structure calculations have been reported for cubic and orthorhombic SrRu03tJ and 
for Sr2RuO40a, using the general potential linearized augmented plane wave method (LAPW). An analysis of the 
resulting electronic structure was given in Refs. [|J|. At the moment, essentially all observable properties which could 
be reliably calculated from the electronic structure results agree well with the experiment. On the other hand, there 
are several properties of this material that cannot be reliably obtained directly with the LAPW method due to the 
high computational cost. At the same time, there is substantial interest among both theorists (as a starting point for 
many body models) and experimentalists in having_a.simple, fast tight binding model Wjhich would accurately describe 
the LDA electronic structure near the Fermi levelrQ. Existing tight binding modelsErH do not reproduce accurately 
enough the details of the fermiology, e.g. the degree of nesting. For this reason we have performed highly accurate 
O ■ tight-binding parametrizations of this system by fitting to LAPW results. 
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For the cubic SrRu03 we constructed orthogonal tight-binding (TB) Hamiltonians by fitting to the spin-polarized 
band structure of Ref. [l]. We followed a procedure similar to the one used for other oxide perovskitesu, based on 
14 LAPW bands for each of 165 k-points in the irreducible Brillouin zone. Our TB Hamiltonian results in a 14x14 
secular equation based on the d-Ru and p-0 orbitals. This Hamiltonian contains 11 first and 2 second nearest neighbor 
Slater-Koster parameters (listed in Table |) for each spin determined by least-squares fitting to the LAPW results. 
The overall rms error for all 14 bands was 14 mRy. This high level of accuracy is demonstrated in our energy band 
and densities of states figures discussed in the next section. As can be seen from Table |, we have included the crystal 
field splitting e g — t^g on the Ru d-states, which is 3.34 eV (T) and 3.23 eV(j). There is also significant crystal field 
splitting — p a on the oxygen site, 1.33 eV (f) and 1.41 eV(|) (a similar effect exists in manganite perovskites, cf. 
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Rcf. [l(]). The exchange splitting for Ru is relatively small: 0.57 eV fag) and 0.45 eV (e g ). We also have exchange 
splitting on the O sites (0.12 eV for p a and 0.20 eV for p n ), in accord with the general role of Ru-0 hybridization 
in ruthenates, as discussed in Ref. ||. The strongest interactions in the system are (Table |) the Ru-0 pda and pdir, 
as expected. We also note that in order to obtain a good fit we had to include the second nearest ppa and ppir 
interactions. Another observation which can be made from Table || is that the rigid-band Stoner picture describes the 
physics of magnetic splitting rather well: there is hardly any difference between the hopping intergal in the spin up 
and in the spin down channels. 

For Sr2RuC>4, we followed a similar procedure to construct the orthogonal Hamiltonian. The details are as follows. 
We fitted to the LAPW results of Ref. | using 18 bands for each of 207 k-points in the irreducible Brillouin zone. 
The size of our Hamiltonian matrix became 27x27 based on the d— Sr, d— Ru and p— O orbitals. Having in mind 
applications to superconductivity theory, we aimed at having a virtually perfect fit in the window of ±1 eV around 
the Fermi level. For that inclusion of the d— Sr orbitals was unavoidable. This Hamiltonian contains 30 first and 2 
second nearest neighbor parameters. Since we used an orthogonal TB model, the physical nonorthogonality between 
d-Ru and p-0 orbitals (which have bigger overlaps than other orbitals in the system) spills over into nonzero second 
nearest neighbor parameters in an orthogonal Hamiltonian, specifically Ru-02 pda, and 01-02 ppa, as well as into 
direct hopping integrals between Ru (the fact that the fitted values for Ru-Ru hopping do not come from actual 
overlap of the wave functions is corroborated by the fact that the ddS parameters have a far stronger effect on the 
band structure than the dda and ddir ones). In our notation 01 and 02 describe the in-plane and apical oxygens, 
respectively. The parameters are listed in Table H. The overall r.m.s. error for this fit is « 10 mRy for all 18 bands 
and is less than 2 mRy for the bands crossing the Fermi level (13th to 15th), for which we used a higher weight in 
order to reproduce the Fermi surface very accurately. 



RESULTS 

In Fig. [l] we show the TB energy bands of cubic SrRu03 for both spin up and spin down. We note that the TB 
bands, especially near the Fermi level (Ep) are almost indistinguishable from the LAPW ones. In the (100) direction 
we note near Ep a flat band that is occupied for spin- up and above ~Ep for spin-down. We also see a small electron 
pocket for spin-up and a gap for spin-down centered at T. We observe a hole pocket for spin-down centered at M that 
vanishes for spin-up. In Fig. || we show the TB DOS for both spins. In agreement with the LAPW results, we find a 
pronounced peak just below Ep for spin-up and just above Ep for spin-down, while DOS at Ep is approximately the 
same for both spins. Both these peaks are a mixture of t2 g -Ru and p^-0 states as can be seen from the decomposition 
of the DOS, in accord with the discussion in Ref. [|. On the other hand, the e g contribution to DOS at the Fermi 
level is very small for spin up and nearly zero for the spin down channel. 

In Fig. |3| we show a comparison of TB and LAPW energy bands for Sr2Ru04. The bands above 2 eV were 
not fitted. In the low bands far away from Ep there are minor discrepancies. However, the bands near Ep display 
excellent agreement between TB and LAPW. This also results in a virtually perfect match of the Fermi surfaces (Fig. 
|j). The degree of the nesting of the pseudosquare Fermi surface sheets is reproduced very well, which is crucial for 
the calculations of susceptibilities and response functions (cf. Ref. |ll|). Also the (very weak) k z dispersion, which 
manifests itself through small violations of the mirror symmetry in Fig. ||, is reproduced. 

Apart from computational efficiency, TB models have another advantage with respect to the first principles methods: 
electronic velocities, which are usually calculated in the first principles methods by numerical differentiation of electron 
eigenenergies, can be obtained analytically from the TB Hamiltonian without any loss of accuracy: 



dE kQ I 

Vk « = w = \ ka 

R 



dH 



(1) 



where summation is performed over all neighbors for which hopping integrals were included in the Hamiltonian. The 
first line follows from the Hellman-Feynman theorem, and the second from the TB expression for H = X)R e * kR ^R- 
The advantage of using this formula is particularly clear when there are many band crossings in the energy range 
of interest, in which case numerical differentiation requires an extremely fine mesh. This is even more true for the 
physical quantities that depend on an electronic mass (second derivative of the eigenenergy) , like Hall coefficient. 
Since the error of numerical differentiation increases rapidly with the order of the calculated derivative, eliminating 
the second differentiation by using the Hellman-Feynman velocities is particularly helpful in this case. To illustrate 
this approach, we computed a quantity that defines the vortex lattice geometry in superconducting Sr2RuO40 and 



2 



depends in a non-trivial way on the electronic velocities. AgterbergHI showed that this geometry is directly related to 
the symmetry of the order parameter. For an isotropic triplet p-wave pairing, one expects a triangular vortex lattice. 
For a d wave pairing a square lattice appears, oriented either along (100) or along (lljQ) directions. It was observed 
experimentallyE] that in Sr2Ru04 Abrikosov vortices form a square lattice; AgterbergLl performed calculations for a 
model of a 2D triplet superconductor, where the angular anisotropy comes entirely from the fermiology. He derived 
a criterion that controls the geometry of the vortex lattice: 



2{5v*f-{v*f _(vl-v, 



»= ^ ' = 2 I 2 -1- (2) 



2 +v 2\ 



The absolute value of this quantity determines whether the vortex lattice will be triangular or square (for v > 0.014 
it is always square), while its sign defines the orientation of the lattice when it is square. If v < it is aligned with 
the crystal lattice, as it is observed in the experiment. 

Agterberg pointed out that when order parameters are different for the three bands crossing the Fermi level, the 
actual vortex lattice will be defined by the sign of the weighted v parameter, which in turn depends on the order 
parameters A in all three bands, 

„„, = (y a (v*)l Al + Vf) {v% Al + v, (vX ^)li(v 2 )l K + (v% A 2 + A*). (3) 

He estimated, using a 3-band TB model of Ref. ||, that v a w Vp w 0.6, while rs —0.6. This calculation is to be 
taken with a grain of salt, though, as the model used for the band structure is rather approximate. However, if indeed 
v a ~ vp ~ — ^7 > 0, even a modest interband anisotropy of the order parameter such as that computed in Refs. |^JTl| 
could ensure agreement with the experiment. Of course, since Eqs. EL[| involve complicated averages over the Fermi 
surface, the result is very sensitive to the accuracy of the Brillouin zone integration. 

In Fig.| we show the results of our calculation of v, using the present TB model, Eq.|l], and the tetrahedron 
integration technique with the 1507 irreducible k-points in the Brillouin zone. First, we observe that even with such 
a fine mesh and relatively simple Fermi surface the numerical differentiation of eigenenergies produces noisy curves, 
compared to the analytical differentiation. Second, we find that there is little dependence of v on the Fermi energy 
in the range ±150 meV, although the van Hove singularity is located only 70 meV above the Fermi energy. Third, 
we see that the calculated values of v (v a — 0.42, va = 0.70, v n = —0.45) are not too far from Agterberg's estimates. 
Fourth, we can use Eq. || to estimate the lower limit on the ratio A 7 /A Q , assuming A a = Ap. In order to reproduce 
the experimental observation that v e ff < 0, and using our calculated ratios (v 2 ) a : (v 2 ) : (v 2 ) = 1.0 : 0.55 : 1.0, 
we find that A 7 /A Q > 1.2, which is to be compared, for instance, with A 7 /A Q «2.2 in Ref. 

Another interesting Fermi surface property, which requires highly accurate integration over the Brillouin zone, 
is the Hall coefficient. It depends on the second derivatives of the eigenvalues (see Appendix), which makes nu- 



merical different iati on undesirable. We used our -IB fit to calculate the Hall coefficient both in diffusive (Eq.A.3) 



and ballistic(Eq. A4) regimes. Experiments showed that the Hall coefficient is negative (electron-like) at T — » 0, 
R = —1.15 x 10~ lu m 3 /C, but grows rapidly and changes sign at T » 35 K, and stays close to zero at the higher 
temperatures. It was suggested that this is a manifestation of the different temperature dependence of the relaxat ion 



time in different band Alternatively, one may think that at T ~ 35 K a crossover between the Eq.b\4| and Eq.|A3| 
takes place. 

We checked these hypotheses by calculating the conductivities a xy and a xx (see Appendix) for all three Fermi 
surface sheets. In both cases they are calculated up to a common constant factor (constant relaxation time r or 
constant mean free path A). The results were, in atomic units, as given in Table E9. First of all, we observe that the 
difference between the two approximations is very small. Second, we see that the net Hall conductivity originating 
from the xz/yz sheets of the Fermi surface (a and /?, in standard notations) is electron-like. Thus the Mackenzie et 
al explanation of the sign change at T ~ 35 K requires a substantial difference of the temperature dependence of r a 
and T/3 (not only of T a ^p and r 7 ), which is physically hard to support, considering the origin of these bands. 

However, as noted in a previous papera, to fit the experimental de Haas-van Alphcn extremal area exactly, one has 
to shift the LDA bands by a few mRy up (a) or dowp-.(/3, 7). In order to investigate this, we applied shifts of 5, -4, 
and -1.5 mRy to the bands a, (3, and 7, respectivelytil This procedure reduces of the energy distance between the 
van Hove singularity and the Fermi level from w 60 to « 20 meV, which makes the difference between the ballistic 
and the diffusive a xy more pronounced. The results are shown in Tabic IV. One immediately observes that now 



the net Hall conductivity in the bands a and (3 is positive. Thus if the mean free path in those two bands changes 
with temperature slower than that for the 7 band (to some extent, in spirit of the "orbital-dependent" conjecture of 
Agterberg et acB), one may, in principle, expect a sign change of the net Hall coefficient. In fact, one needs to assume 
that the ratio Aq,,^/A 7 increases from w 1 at T = to ks 4 at T ss 40 K to explain the observed sign change. This 
would be unusual, but certainly not impossible. 
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CONCLUSIONS 



To conclude, we presented highly accurate tight binding fits for cubic SrRu03 and tetragonal Sr2Ru04 perovskites. 
Using nearest neighbor and a few selected next nearest neighbor hopping parameters we were able to reproduce 
the first-principles band structure of Sr2Ru04 in the physically relevant energy range near the Fermi level with the 
accuracy of a few mRy This allows for fast generation of the electronic eigenenergies at a superfine mesh in the 
Brillouin zone, and fast and virtually exact calculations of the electronic velocities at the same mesh. This powerful 
technique provides a possibility of investigating physical properties with high sensitivity to the details of the band 
structure and Fermi surface topology. 

We report two application of the present fit to actual physical problems: First, we calculated the anisotropy 
parameter v relevant for the superconducting vortex lattice geometry, and derived from that the lower bound for the 
interband anisotropy of the order parameter: to ensure compatibility with the experiment, the order parameter in 
the xy band should be at least 20% larger than in the yz/zx bands. Second, we calculated Hall conductivity for the 
three bands crossing the Fermi level, both in the constant-relaxation-time and in the constant-mean-free-path regimes. 
Having these numbers we tested the two most natural hypotheses for the temperature-induced sign reversal of the 
Hall conductivity: crossover between the two above-mentioned regimes and different temperature dependence of the 
xy— and yz/zx— carrier mobility. We can confidently rule out the former hypothesis. The latter hypothesis cannot 
be definitely ruled out; for the Fermi surface adjusted to fit the de Haas - van Alphen cross-section the desired sign 
reversal can be obtained in the constant-mean-free-path regimes, if the mean free path for the xy— carriers changes 
between O K and 40 K at least four times faster than that for the yz/zx— carriers. A possible reason for such a 
behavior would be substantially stronger electron-paramagnon coupling for the xy— electrons. 
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APPENDIX: DERIVATION OF THE HALL FORMULAS 



It is convenient to use the 2D expression given by OngB 

e 2 B 



d 
dk„ 



v x (k)T 



ki 



ki' 



(Al) 
(A2) 



ki 



which defines the Hall coefficient, R = —a xy /Ba 2 x . Here VI is the unit cell volume, and summation is over all bands 
and all states in the Brillouin zone. In the constant-relaxation-time approximation this reduces to the following 
expression: 
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Rd = — ^ 5(Eki - E F )(vlHyy ~ V x Vy[l xy )/2 



Y,6(Eki-E F )v 2 x 



(A3) 



where /i xy (k) = d 2 E^/ 'dk x dk y etc. As noticed by Mackenzie et a0, at very low temperature a more appropriate 
approximation is where the mean free path, Tk is a constant, in which case we have instead 

2 



R b = ^25(E ki - E F )(vlflyy - V x Vy[l X y)V ^ / 2 



ki 



Y,S(E kl -E F )v x 



ki 



(A4) 



1 D.J. Singh, J. Appl. Phys. 79, 4818 (1996). 



4 



2 D.J. Singh, Phys. Rev. B 52, 1358 (1995). 

3 T. Oguchi, Phys. Rev.B 51, 1385 (1995). 

4 I.I. Mazin and D. J. Singh, Phys. Rev. B56, 2556 (1997). 

5 LI. Mazin and D. J. Singh, Phys. Rev. Lett. 79, 733 (1997). 

6 A.I. Liechtenstein, L. Liebsch, unpublished. 

7 D.F. Agterberg, Phys. Rev. Lett. 80, 5184 (1999); Phys. Rev. B58, 14484 (1998). 

8 C. Noce and M. Cuoco, Phys. Rev. B59, 2659 (1999). 

9 J. P. Julien, D.A. Papaconstantopoulos, D.J. Singh, W.E. Pickett, F. Cyrot-Lackmann, Physica C, 220, 359 (1994). 

10 D. A. Papaconstantopoulos and W. E. Pickett, Phys. Rev. 57, 12751 (1998). 

11 I.I. Mazin and D. J. Singh, Phys. Rev. Lett. 82, 4324 (1999). 

12 T.M. Riseman et al. Nature 396, 242 (1998). 

13 A. P. Mackenzie et al, Phys. Rev. B54, 7425 (1996). 

14 The last number differs from the one in Ref. ^ (3 mRy) because of the difference between the TB fit and the LAPW results. 

15 D.F. Agterberg, T.M. Rice, and M. Sigrist, Phys. Rev. Lett. 78, 3374 (1997). 

16 N.P. Ong, Phys. Rev. B43, 193 (1991). 



5 




FIG. 1. TB band structure of cubic SrRuOa for two spin directions. 
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FIG. 2. TB DOS of cubic SrRuQ 3 for the two spin directions. 
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FIG. 3. Band structure of Sr 2 Ru0 4 in LAPW (circles) and in the TB fit (lines). 
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FIG. 4. Fermi surface of S^RuCU in LAPW (solid lines) and in the TB fit (dashed lines) 
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FIG. 5. Anisotropy parameter v (see text) as the function of the distance to the Fermi level. Solid lines are results of the 
integration with exact TB velocities, broken lines use velocities from tetrahedron linear interpolation. In both cases 30x30 x 10 
mesh was used. 
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TABLE I. Slater-Koster parameters for the cubic SrRuOa (mRy). 



spin 


Ru onsite 

tig e g 


O onsite 

Pa P-k 


Ru-Ru 
dda ddn ddS 


0-0 

ppa ppn 


Ru-O 
pda pdn 


0-0 

ppa a ppn a 


T 
1 


493 738 
534 772 


307 405 
316 420 


-28 -12 1 
-27 -11 1 


43 -11 
43 -11 


168 -94 
165 -97 


14 -4 
14 -5 



a second nearest neighbors 



TABLE II. Slater-Koster parameters for Sr2Ru04 (mRy). 



Ru onsite 


01 onsite 


02 onsite 


Sr onsite 


Ru-Ru 


Ol-Ol 


02-02 


01-02 


xy yz, zx e g 
487 516 659 


P 

338 


z x,y 
392 436 


xy yz, zx e g 
1032 1318 1874 


dda ddn ddS 
-7 -3 -6 


ppa ppn 
37 -7 


ppa ppn 
7 4 


ppa ppn ppa a 
44 -18 -7 



Sr-Sr 


Sr-Ru 


Ol-Sr 


Ol-Ru 


02-Sr 


02-Ru 


ddn 
39 


dda 
-80 


ddn 
-7 


ddS 
-38 


pda pdn 
-47 110 


pda pdn 
190 96 


pda pdn 
61 -58 


pda 
160 


pdn 
68 


pda a 
17 



"second nearest neighbors 



TABLE III. Conductivities and Hall coefficients in atomic units. The volume of the unit cell Q. =95.4 A 3 , r and A are the 
relaxation time and the mean free path, respectively. Atomic units for R are such that for a parabolic band R — 1/n, where n 
as given in the first columns is the total number of electrons (a, 7) or holes (/3) in the corresponding bands. 



band 


n 


4JtQ 






a b x J\U 


R d O. 




Rcx P Q(T — 


0) 


Q 


-1.02 


2.0 


6.0 


-4.0 


-36 


-0.95 


-0.98 






P 


0.28 


1.1 


3.2 


3.8 


32 


3.2 


3.2 






7 


-1.26 


2.0 


6.6 


-2.8 


-30 


-0.65 


-0.68 






total 


-2.00 


5.1 


15.2 


-3.0 


-34 


-0.12 


-0.15 


-0.20 






TABLE IV. 


The same as Table 


|lV|, but for the bands shifted to match de Haas 


van 


Alphen cross section areas. 




band 


n 


ct xx /tQ, 


a b xx /\Q, 


aiy/rQ. 


a xy /\Q, 


R d n 


R b n 


Rcxp&(T — 


0) 


a 3 


-0.92 


2.0 


5.7 


-4.0 


-35 


-1.03 


-1.05 






(31 


0.25 


1.0 


3.0 


3.4 


37 


3.4 


4.17 






72 


-1.35 


1.9 


6.7 


-1.8 


-29 


-0.52 


-0.64 






total 


-2.02 


4.9 


15.4 


-3.4 


-31 


-0.14 


-0.13 


-.20 
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